The Writhe Distribution Of Stretched Polymers 
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Motivated by experiments in which single DNA molecules are stretched and twisted we consider 
a perturbative approach around very high forces, where we determine the writhe distribution in 
a simple, analytically tractable model. Our results are in agreement with recent simulations and 
experiments. 
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In recent years, there has been much interest in the 
area of statistical mechanics of semiflexible polymers. 
These studies have been motivated by experiments [1] 
on biopolymers in which single molecules are stretched 
and twisted to measure elastic properties. These experi- 
ments are designed to understand the role of semiflexible 
polymer elasticity in, for instance, the packaging of these 
polymers in a cell nucleus. Twist elasticity plays an im- 
portant role in several biological functions. The first step 
in packaging DNA in a cell nucleus a few microns across 
involves DNA-histone association which makes use of su- 
percoiling in an essential way. The process of DNA tran- 
scription can generate and be regulated by supercoiling 
[2]. Here we focus on a particular class of experiments 
which probes the twist elasticity of DNA. 

In the experiments of Strick et al [2] the ends of a single 
molecule of double stranded DNA are attached to a glass 
plate and a magnetic bead. Magnetic fields are used to 
rotate the bead and magnetic field gradients to apply 
forces on the bead. By such techniques the molecule is 
stretched and twisted and the extension of the molecule 
is monitored by the location of the bead. One thus gets 
the extension of the molecule as a result of the applied 
twist and force. 

The simplest theoretical model used to interpret the 
experiments makes use of the fact that the molecule is 
under such high tension that it is essentially straight. 
We call this limit the paraxial limit of the elasticity of a 
semiflexible polymer keeping the optical analogy in mind 
[3,4]. In this situation, the molecule, being straight and 
taut, cannot intersect itself and so one does not expect 
self- avoidance effects to be important. In computing the 
partition function one simply sums over all configurations 
without regard to self-avoidance. This model is instruc- 
tive because it is analytically tractable and enables us for 
the first time to derive a simple explicit analytic expres- 
sion for the writhe distribution. This is the central result 
of our analysis. 

Before discussing the schematics of the derivation of 
the writhe distribution let us define a few pertinent quan- 
tities: link, twist and writhe. For an open polymer of 



the kind that is used in the twist-stretch experiments 
one can define link {Lk), an arbitrary real number, as 
the externally imposed twist given by 27rn where n is the 
number of applied turns on the bead. Twist {Tw) corre- 
sponds to the integrated rotation of the polymer around 
its backbone and writhe [Wr) pertains to the twist of the 
polymer backbone and is captured by the rotation of the 
tangent vector. In computing the writhe distribution we 
make use of the fact that link Lk is related to twist and 
writhe via Lk — Tw + Wr [5-7]. It is convenient to go 
to the conjugate space and work with the variable B, the 
generator of link Lk. In this space the partition function 
Z{B,f) neatly factors into writhe Zw{B,f) and twist 
ZxiB,/) [8]. The distribution Z{Lk,f) is a convolution 
of the writhe distribution and the twist distribution. We 
compute the link distribution in the a = Lbp/Ltp (the 
ratio of the bend and the twist persistence lengths) 
limit. In this limit the writhe distribution coincides with 
the link distribution. This is simply because in the a ^ 
limit it is very expensive to twist the polymer around a 
straight backbone and the polymer goes into a bending 
mode resulting in a twisting of the polymer backbone 
[5] . From the writhe distribution computed in this man- 
ner, one can recover the link distribution by convolving 
it with the distribution of "dynamical twist" . 

Our starting point is the Worm Like Chain [WLC) 
model in which the polymer is modelled as a framed space 
curve C — {x{s), ei{s)}, i — 1,2, 3, where < s < L is the 
arclength parameter along the curve. The unit tangent 
vector 63 = dx/ds to the curve describes the bending 
of the polymer while the twisting is captured by a unit 
vector ei normal to 63. In keeping with the optics analogy 
[3,4], we refer to ei as the polarisation vector. 62 is then 
fixed by 62 = 63 x ei to complete the right handed moving 
frame ei{s), i = 1,2, 3. The energy £[C\ of a configuration 
of the polymer is a sum of contributions coming from its 
bending and twisting modes. 

In the presence of large forces |i^| — > 00, the molecule 
is stretched taut and there is an energy barrier for the 
polymer to pass through itself. In fact, the molecule is 
constrained to lie in nearly a straight line between its 
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ends with small deviations. As mentioned earlier, in 
this regime the polymer being essentially straight cannot 
cross itself and thus self-avoidance effects present in a real 
polymer are automatically taken into consideration. Un- 
der these conditions the tangent vector only makes small 
deviations from the z direction. We can approximate the 
sphere of directions [5] by a tangent plane at the north 
pole of the sphere. We call this limiting model of the 
WLC the paraxial worm like chain or the PWLC model. 

In this straight taut limit the polymer Hamiltonian 
[5,8-10] reduces to : 

HpwLC - y + 2^ + — - /(I - y) 



= Hp-f + 



B^a 



where Hp is the Hamiltonian of interest in the parax- 
ial limit after we take out a constant piece, a is the 
ratio of the bend persistence length Lbp and the twist 
persistence length Ltp- The constant B corresponds to 
the conserved momentum conjugate to the Euler angle 
V'. / = FLbp/UbT where F is the stretching force 
and kpT is the thermal energy. The 'vector potential' 
A^ = B^. Thus the PWLC maps on to the problem of 
a particle moving on a plane in the presence of a magnetic 
field B and an oscillator confining potential which arises 
from making a small expansion for the stretching force 
-/cos^(-/cos^«-/(l-f ) = -/ + /f ) [8,10]. No- 
ticc that, in contrast to the regime of low tension [5], in 
this high tension regime the polymer cannot release an 
imposed twist by passing through itself because of the 
condition of suppression of configurations in which the 
polymer folds back onto itself. This implies that in con- 
trast to the WLC model [5], in the PWLC model the free 
energy, torque-twist relation and other related distribu- 
tions are not periodic functions of the imposed twist. 

Introducing Cartesian coordinates = 9cos(p and 
^2 = sin on the tangent plane at the north pole 
of the sphere of directions one can express the small 9 
Hamiltonian Hp as follows: 

Hp = - + - A^^ f + ^ ^^2^ 

(1) 

where A^^ = —5^2/2, = The corresponding 

partition function Z can be written in terms of the eigen- 
values En and eigenfunctions of HpwLC as follows: 



Z{B,f,(o,{L)=Y.' 



<(Co)Wn(^L) 



where = (6(0), 6(0)) and 6 = (6(i),6(i)) are 
the initial and final tangent vectors at the two ends 



of the polymer. In order to simplify our analysis fur- 
ther we confine ourselves to the limit of very long poly- 
mers. Many of the experiments involving biopolymers 
such as DNA explore this limit of very long polymers 
which is also theoretically more tractable. For long 
polymers (/? = — > oo) only the lowest eigenvalue 

EoiBJ) = ^f + ByA-f+^ [10,11] oiHpwLC dom- 
inates the expression for the partition function. Thus the 
partition function can be written as: 



ZiBJ)0^, 



-0Eo{BJ) 



(2) 



Writhe Distribution: 

Let us consider the link distribution Z{Lk, /): 



(3) 

where the phase (/'(B, Lk, /) is given by 4>{B, Lk, /) = 
iPEo{B, /) — jSB^^. Since we are working in the limit of 
long polymers we can compute this partition function us- 
ing the stationary phase or saddle point method [8] where 
only the stationary value (/>(i?sf , LA;, /) of the phase dom- 
inates. <l){Bst, Lk, f) is the central quantity of interest 
from which all the relevant elastic properties of the taut 
polymer can be derived. A similar perturbative analysis 
was done by Moroz and Nelson [10], who in fact carry 
the analysis out to higher orders in perturbation theory. 
What is new in our treatment is an explicit analytical ex- 
pression for the writhe distribution in the straight taut 
limit. 

Here we outline the derivation for the writhe distribu- 
tion. The condition for stationarity (^ = 0) satisfied 
by the stationary value B^t of B is: 



[ilk - aBst) 



(4) 



where Lk = We restrict to the case of a = 0. In 
this case the equation simplifies and we get the following 
stationary value of Bgt- 



±i 



[1 + 4lV] 



Of the two roots, only the positive root is the physically 
relevant one consistent with the saddle point approxima- 
tion. Setting iBst = t. where r has the interpretation of 
torque, we get the following Torque-Link relation: 



T = 



^J[1 + ALk^] 



(5) 



Notice that for small Lk we get a linear Torque-Link 
(r — Lk) relation which goes over to a Torque-Link re- 
lation independent of Lk in the limit of large Lk. This 
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is consistent with recent experiments [12] and mimeri- 
cally generated plots [8]. Inserting the expression for 
the stationary value B^t of B into the partition func- 
tion Z{Lk, /) we get the pertinent writhe distribution 
P{W, /). To compute this distribution we have made use 
of the fact that for a — > 0, twist is extremely expensive 
and the applied twist goes completely into the bending 
mode. Thus, the link distribution in this limit corre- 
sponds to the distribution P{W, /) of writhe. Given the 
writhe distribution P(W, /) obtained in this manner the 
link distribution P{Lk, f) can be constructed for all val- 
ues of a by convolving it with the twist distribution. In 
the generating function space one simply needs to multi- 
ply the writhe partition function Z\Y{B,f) by a simple 

Gaussian factor Zt{B, /) — e 2 pertaining to the pure 
twist distribution at finite a. 

The analytic form of the scaled distribution P{W) = 
^P(fi 'f ) writhe and stretch in the high tension regime 
(See Fig. 1 ) is: 



Piw) = exp[-py/f{vT+m^ - 1}] 



(6) 



Here / = with F the applied stretching force. This 

analytic form of the writhe distribution is the central 

result of this paper. Plots of this distribution are dis- 
played in Fig. 1. The form reduces to a Gaussian form 
{P{W) w exp[-2/3v/7W^]) for small values of writhe W 
and goes over to P{W) « exp[— 2/3^/71^^1] for very large 
values of writhe W. The writhe distribution (Eq. (6)) 
represented in Fig. 1 has all the expected features- it 
peaks near smaller values and dies off for larger values of 
writhe. Writhe gets suppressed with increasing strength 
of the stretching force. These qualitative features are 
in agreement with recent simulations [13] of writhe as 
a function of the stretching force. The explicit expres- 
sion for the writhe distribution presented here is exact in 
the high tension limit and we expect quantitative agree- 
ment between the predicted distribution and future ex- 
periments probing the writhe distribution in this regime. 

To summarize, we have, for the first time obtained an 
explicit analytical expression for the writhe distribution 
of a semiflexible polymer in the high tension regime. The 
expression for the writhe distribution is simple and trans- 
parent and the qualitative features agree well with avail- 
able computer simulations [13]. For very large forces a 
DNA molecule undergoes force induced denaturation [2] 
and therefore the distribution predicted here may not be 
directly applicable to DNA experiments at very high ten- 
sion. However, one can test the predictions against exper- 
iments with other semiflexible polymers. We therefore, 
expect this work to generate interest amongst experimen- 
talists to measure the writhe distribution of stretched 
polymers. The distribution computed here is also rele- 
vant to depolarized light scattering in turbid media in 
the limit of small angle scattering [3,4,14]. We also have 
an explicit analytic form for the Torque- Twist relation 



which is in agreement with recent experimental data [12]. 
In this model the mean-squared writhe fluctuation which 
corresponds to the second derivative of the conjugate 
distribution Z{Lk, f), diverges at t = 2^//, the point 
at which a 'buckling instability' sets in for the polymer. 
More explicitly, the divergence of the second moment of 
the writhe distribution at the buckling instability point 
T = 2y/f has the following form which can be tested 
against future experiments: 



< >= 



f 



4[/-^]3/2 



The second moment of the distribution is consistent with 

earlier predictions for the mean squared writhing angle 
for long tense molecules [15]. At the buckling instability 
there is a divergence of the writhe fluctuations. Since 
such a divergence makes the polymer backbone fluctu- 
ate violently one expects it to lead to a corresponding 
divergence in the mean squared extensional fluctuations 
< >. This has, in fact been probed in some recent 
experiments [16]. In the 'paraxial' limit we find that 

<W^> , /■ . 
The 'paraxial' approximation breaks down for large 

values of the applied twist (more precisely for r > 
[See Eq. 4] ) in which case the polymer explores configu- 
rations which deviate considerably from the straight taut 
limit. In future we would like to explore the low tension 
nonlinear regime for the writhe distribution where phe- 
nomena like plectoneme formation would play an impor- 
tant role and nontrivial self-avoidance effects [3,15,17,18] 
need to be taken into consideration. The present work 
will provide a limiting check on calculations done in the 
nonlinear regime. As we mentioned earlier, the writhe 
distribution has important implications in the context of 
transcription and gene regulation. Therefore, a complete 
understanding of writhing of a biopolymer backbone and 
its stabilization is of relevance to current research. 




FIG. 1. The Writhe distribution P{W) for / = 2 and 
/ = 5(dashed curve) for = 10. 
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